ON THE DISTRIBUTION OF CARMICHAEL NUMBERS 



ARAN NAYEBI 

Abstract. Erdos conjectured in 1956 that there are x^~°^^^ Carmichael numbers up to x. Pomerance 



1- 



{l + o(l)} logloglogx 



made this conjecture more precise and proposed that there are x log log ^ Carmichael numbers 

up to X. At the time, his data tables up to 25 • 10^ appeared to support his conjecture. However, Pinch 
extended this data and showed that up to lO'^^ , Pomerance's conjecture did not appear well-supported. Thus, 
the purpose of this paper is two-fold. First, we build upon the work of Pomerance and others to present 
an alternate conjecture regarding the distribution of Carmichael numbers that fits proven bounds and is 
better supported by Pinch's new data. Second, we provide another conjecture concerning the distribution 
of Carmichael numbers that significantly sharpens Pomerance's heuristic arguments. We also extend and 
update counts pertaining to pseudoprimes and Carmichael numbers, and discuss the distribution of One- 
Parameter Quadratic-Base Test pseudoprimes. 



1. Introduction 

Fermat's "little" theorem states that if b is an integer prime to n, and if n is prime, then 
(1.0.1) 6" = & (mod n). 

When gcd(6,7i) ~ 1, we can divide by b, 

(1.0.2) 6"-i = 1 (mod n). 

A composite natural number n for which 6"^^ = 1 (mod n) for any fixed integer 6 > 2 is a base b pseudoprime. 
A positive composite integer n is a Carmichael number if = 1 (mod n) for all integers b > 2 with 
gcd(6, n) = 1. The importance of Carmichael numbers is that they test the limits of the Fermat primality 
test, which ultimately led mathematicians to formulate more effective tests. Furthermore, there is little that 
is known about them; for instance, the infinitude of Carmichael numbers has only recently been proven by 
Alford, Granville, and Pomerance [15] . 

Let ^b{x) denote the number of base b pseudoprimes < x and let C{x) denote the number of Carmichael 
numbers < x. In 1899, Korselt [20] provided a method for identifying Carmichael numbers 

Theorem 1.0.1. An odd number n is a Carmichael number iff n is squarefree and p — 1 | n — 1 for all p \ n, 

where p is a prime number. 



As a consequence of Theorem 11.0. 11 it is easy to see that Carmichael numbers have at least three prime 
factors. 

In 1910, Carmichael 22J found the smallest Carmichael number to be 561 = 3 • 11 • 17. 

Based on Korselt's criterion, Erdos Tf] formulated a method for constructing Carmichael numbers, which 
was the basis for the proof of Alford, Granville, and Pomerance [15 . His notion was to replace "p—1 \ n — 1 
for all p I n" in Theorem 1 1 . . 1 1 with L | n — 1 for L lcmp|„(p — 1). By focusing primarily on L, Erdos 
found every p for which p — 1 \ L and then tried to find a product of those primes in which = 1 (mod L) [8\ . 
His results heuristically suggested that for sufficiently large x, 

(1.0.3) C(x) =a;i-°(i). 

More convincingly. Theorem 4 of [TS] shows that (|1.0.3p holds if one assumes widely-believed assumptions 
regarding primes in arithmetic progressions. However, drawing upon available data at the time. Shanks [23] 



Date: February 2009. 

2000 Mathematics Subject Classification. Primary 11Y35, 11N60, Secondary 11N05, 11N37, 11N25, llYll. 

1 



2 



ARAN NAYEBI 



was skeptical of (jl.0.3p because the counts of Carmichael numbers seemed to have noticeably fewer prime 
factors than those predicted by Erdos' heuristic. 

Granville and Pomerance [8 conjectured that the reason for the difference between the computational 
evidence and the argument of (|1.0.3p stems from a grouping of Carmichael numbers into two distinct classes, 
namely primitive and imprimitive. If we let g = g{n) := gcd(j»i — l,p2 ~ 1, • ■ ■ iPfc ~ 1) for a squarefree 
integer n = piP2 ■ ■ 'Pk and put pui = pt — \ for some integer aj, then n is a primitive Carmichael number if 
g{n) < [ai, • • • , a^], and imprimitive if otherwise. Thus, since the observations of Shanks are more applicable 
to imprimitive Carmichael numbers and those of Erdos are more applicable to primitive Carmichael numbers, 
and most Carmichael numbers are in fact primitive whereas most Carmichael numbers with a fixed number 
of prime factors are imprimitive, then the two conjecturers easily reached different conclusions. 

Interestingly, Pinch's counts of fc-prime Carmichael numbers up to 10^^ |16| reproduced in Table [1] imply 
that the number of prime factors of primitive Carmichael numbers tends to increase as x gets larger. As can 

Table 1. Counts of fc-prime Carmichael numbers 



Bound 


Cz{x) 








Cjix) 


Cfi{x) 


C^{x) 


Cio{x) 


Cii{x) 


Cl2[x) 


C[x) 


10^ 


1 





























1 


10^ 


7 





























7 


10^ 


12 


4 


























16 


10^ 


23 


19 


1 























43 


10^ 


47 


55 


3 























105 


108 


84 


144 


27 























255 


109 


172 


314 


146 


14 




















646 


lOlo 


335 


619 


492 


99 


2 

















1547 


10" 


590 


1179 


1336 


459 


41 

















3605 


1012 


1000 


2102 


3156 


1714 


262 


7 














8241 


1013 


1858 


3639 


7082 


5270 


1340 


89 


1 











19279 


1014 


3284 


6042 


14938 


14401 


5359 


655 


27 











44706 


1015 


6083 


9938 


29282 


36907 


19210 


3622 


170 











105212 


1016 


10816 


16202 


55012 


86696 


60150 


16348 


1436 


23 








246683 


lOi^ 


19539 


25758 


100707 


194306 


172234 


63635 


8835 


240 


1 





585355 


1018 


35586 


40685 


178063 


414660 


460553 


223997 


44993 


3058 


49 





1401664 


1019 


65309 


63343 


306310 


849564 


1159167 


720406 


196391 


20738 


576 


2 


3381806 


1020 


120625 


98253 


514381 


1681744 


2774702 


2148017 


762963 


114232 


5804 


56 


8220777 


1021 


224763 


151566 


846627 


3230120 


6363475 


6015901 


2714473 


547528 


42764 


983 


20138200 



be implied from Table [H for the maximum number of distinct prime factors k{x) <^ iJjf§f^' 
(1.0.4) C(x)^C^(x)+Ci{x)+C^(x) + --- + Cki,^){x), 

where log^-'-' x denotes the j-fold iteration of the natural logarithm for j >2 (we shall use this notation from 
now on). Moreover, if we allow Ck{x) to represent the number of Carmichael numbers < x with precisely 
A; > 3 prime factors, then it is conjectured that 

(1.0.5) Cfe(x) =17fe(a;i/Vlog'x). 

Returning to the Erdos-Shanks controversy, Pomerance [8] sharpened the conjecture in (|1.0.3p for all large 
X in order to be consistent with both Shanks' and Erdos' observations. Define the function h{x) as 

(1.0.6) ^(a;) =x.exp(-/i(x)iH^^^^^|. 

l0g^2J 2; J 

According to Pomerance, distribution of Carmichael numbers is given by 

1 {l + o(l)}log(3) ^ 

(1.0.7) C{x)=x , 
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for X sufficiently large. Unfortunately, according to Pinch 19J, there appears to be no limiting value on h as 
indicated by the recent counts of Carmichael numbers up to lO^^. It is obvious that (|1.0.7p holds iff lim /i = 1 
in (|1.0.6p . However, Pinch [19| explains that the decrease in h is reversed between 10^^ and 10^"', which is 
presented in Table [2] In fact, there is no clear evidence that suggests \m\h = 1. 



Table 2. Values of h{x) 



Bound 


h{x) 


10^ 


2.93319 


10* 


2.19547 


10^ 


2.07632 


10^ 


1.97946 


10^ 


1.93388 


108 


1.90495 


109 


1.87989 


lOi" 


1.86870 


10" 


1.86421 


1012 


1.86377 


10" 


1.86240 


10^4 


1.86293 


10" 


1.86301 


10" 


1.86406 


10" 


1.86472 


10" 


1.86522 


10" 


1.86565 


1020 


1.86598 


1021 


1.86619 



As a result, we present an alternate conjecture 
Conjecture 1.0.2. 

C^{x)^b{x) 



(1.0.8) C{x) 



where !^b.2{x) is the number of two-prime base b pseudoprimes < x and Cz{x) is the number of three-prime 
Carmichael numbers < x. 

From Conjecture 11.0.21 we derive a corollary that the number of Carmichael numbers up to x sufficiently 
large is 



Corollary 1.0.3. 



tp[x2 log X / 



^ dt 



where L{x) ~ exp{ ^ ^ }; "0' = "01 — ' P' '^^'^ ^ '"^'^ primes, and we define 

^a,b,c{p) o,s the number of distinct residues modulo p represented by a,b,c, then the constants C and tt, are 
explicitly given as such, 

(1.0.10) c = 4rV V - ■^)) 

s>l r>s 

gcd(r,s) — 1 
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dim 



2, if4:\m; 



S{m) 



1, i/ otherwise. 

(1.0.11) T3 = K3A, 



K3 



gcd(n, 6) p ' -pf p-a;a.b,c(p) 

n>l p|n a<b<c,n—abc p\n 

p>3 a,b,cpairwise coprime p>3 



S (a, 6, c) 



2, ifa = b = c^O (mod 3); 



1, i/ otherwise. 



Based upon the computation of C made by Galway ,3^ and the evaluation of K3 by Chick and Davies ^3] , 
we beheve that "0' will approach 69.51 and "01 will approach 2.57; although these values are not yet borne out 
by the data. We also demonstrate that CoroUarv ll.O.BI fits the proven upper and lower bounds for C(x), that 
ip' and "01 appear to approach constant values based upon Pinch's data, and we support Conjecture 11.0.21 
through computational efforts. 

In private communication [28j . Pomerance suggests to us that the reason for h{x) not approaching its 
conjectural limit of 1 is that "some secondary terms may be present. So, say in my conjecture, one replaces 
"log^^-* x" with "log'^' X + log^^^cc". It is the same conjecture, since the two are asymptotic... and so the 
Pinch phenomenon is banished". Hence, if secondary terms do indeed exist, then another viable conjecture 
regarding C{x) would be to sharpen the heuristic arguments in [11] which, as a consequence, may better 
match the actual counts of Carmichael numbers. Since these heuristic arguments are dependent upon the 
number of y-smooth numbers up to x, represented by 'i!{x,y), with y in the vicinity of exp{(loga;) 2 }, then 
it would suffice to utilize improvements concerning the asymptotic distribution of these numbers in the 
aforementioned region. As a result of these endeavors, we obtain the more precise heuristic: 

Conjecture 1.0.4. Let tt(x) be the prime counting function, then C{x) for x sufficiently large is asymptot- 
ically close to the function a{x) defined as 

(log(2) a;)2^((logx)i°8' ' exp{-{l + o(l)} log^'^ xlog^^^ x} \ 



(1.0.12) a{x) 



log a; 



The reader may wonder why two conjectures are presented. The reason is that due to the lack of informa- 
tion regarding Carmichael numbers and their distribution, both conjectures are viable to their own merit. 
CoroUarv II. 0.31 asserts that if the values of and ipi can be accurately determined then an asymptotic for- 
mula for C{x) easily follows. Coni ect ure 11.0.41 indicates to us that if secondary terms exist, then the relation 
between the functions 'i'{x,y) and 'i>'{x,y) must be examined further (we explain this concept fully in §3.4) 
to effectively prove an equality for C{x). We should note that the values of C(x) predicted by Corollarv ll.0.31 
and Coniecture ll. 0.41 appear to be closer to the actual values of C{x) than Pomerance's conjecture in (|2.2.2p . 
Moreover, at least up to 10^^, it appears that Conjecture 11.0.41 is presenting more accurate values of C{x) 
than Corollarv ll.0.31 although, this may cease to be the case for larger bounds. 
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Table 3. Comparisons between the actual and predicted Carmichael number counts 



Bound 


L\x) 


69.51x1 






^ {1 + o(1)}1oe(3) ^ 

,-r, 1^-.fT(2) T 
\£ JOgV ' X 


log x-L{x) 


L{x) 


-1 r^^ 
10'' 


1 


301.95 


1092.82 


3.50 


94.89 


10^ 


7 


594.43 


2835.17 


7.81 


365.59 


10^ 


16 


1316.29 


6640.29 


18.18 


1485.33 


10^ 


43 


3131.53 


14806.24 


43.43 


6224.10 


10^ 


105 


7826.17 


32411.27 


107.50 


26636.80 


10* 


255 


20282.91 


71150.56 


274.074 


115803.60 


10^ 


646 


54070.80 


159157.24 


724.86 


509769.35 


10^" 


1547 


147451.71 


367012.00 


1926.56 


2267174.18 


10^^ 


3605 


409716.38 


878601.38 


5245.56 


10171329.99 


1 ni2 
iU 




iiOODo/.OO 


lloobb i .16 


1 /I /I QQ OO 

14400. ZZ 


4oy ( ( u ( y.uy 


1013 


19279 


3309970.24 


5664006.88 


40424.93 


209219668.02 


1014 


44706 


9585268.36 


15162465.67 


114558.014 


957710051.36 


1015 


105212 


28049810.91 


41763706.96 


329251.92 


4407472357.25 


1016 


246683 


82852448.55 


117743387.56 


955940.22 


20382638275.29 


1017 


585355 


246785788.13 


338238941.70 


2796027.81 


94682736406.04 


1018 


1401644 


740679196.52 


986503770.93 


8260103.95 


441642695710.74 


1019 


3381806 


2238429061.23 


2913197684.15 


24637581.64 


2067911761776.64 


1020 


8220777 


6807841639.58 


8692508977.60 


74026750.39 


9717200728399.57 


1021 


20138200 


20826296835.28 


26167265004.43 


224193470.90 


45814162191297.01 



2. Preliminaries 

Before delving into the main results of this paper, we shall first present results regarding pseudoprimes 
and Carmichael numbers that we will explicitly use later on in our derivations. 

2.1. Pseudoprimes. Currently, the tightest bounds for pseudoprime distribution have been proven by 
Pomerance [5] pT| . 

Theorem 2.1.1 (R. A. MoUin 1989, Pomerance 1981). For the base 2 pseudoprime counting Junction, 
exp{(loga;) 207 } < ^2(2^) < x ■ L{x)~ , where L{x) = exp{ ^°^^^^°2) ^ ^ } ■ These bounds are applicable to 
^b{x) for X > xo{b). 

Theorem 2.1.2 (Pomerance 1981). If we allow hin) to denote the exponent with multiplicative order of 2 
modulo n, then n is a pseudoprime (base 2) iff 1 2 (n) \ n — I. 

Conjecture 2.1.3 (Pomerance 1981). The number of solutions w for all n and x sufficiently large is, 

(2.1.1) < X : hiw) = n} < x ■ i(a;)~^+^^''\ lim 9{x) = 0. 

X — >oc 

As a result, the number of base b pseudoprimes for sufficiently large x > xo{b) is conjectured to be, 

(2.1.2) ^tix) ^ X ■ L{xy^ . 

Galway [3] has recently conjectured a formula for the distribution of pseudoprimes with two distinct prime 
factors, p and q, based on a longstanding conjecture of Hardy and Wright concerning the density of prime 
pairs. He noticed that a majority of these pseudoprimes satisfy the relation = where r and s are 
small coprime integers. Thus, we heuristically have 

Conjecture 2.1.4 (Galway 2004). Allow p, q, and d be odd primes, allow ^b,2{x) to represent the counting 
function for odd pseudoprimes with two distinct prime factors, and S^h.iix) := ff\n < x : n = pq,p < 
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q, I^ii(n)}. Hence, as x ^ oo, 
(2.1.3) 
where 
(2.1.4) 

(2.1.5) 
(2.1.6) 

(2.1.7) 



CX2 



log^ X ' 
5{rs)p[rs{r — s)) 



30.03, 



s>l r>s 

gcd(r,s) — 1 



1 -2/d 



J-J- 1 - 1 



(1 - l/df 

d-l 



1.32, 



p(m) = ]J 



2' 



(5(m) = < 



2, ii/4|m; 

1, if otherwise. 



Galway's conjecture is somewhat supported by Table [4] for it appears that C is slowly approaching its 
predicted constant value of 30.03: 



Table 4. Values of C 



Bound 


^6,2(2;) 


C 


10^ 








10^ 


11 


9.331 


10^ 


34 


14.251 


10<^ 


107 


20.423 


10^ 


311 


25.550 


10« 


880 


29.860 


109 


2455 


33.340 


101° 


6501 


34.468 


10" 


17207 


34.908 


10" 


46080 


35.181 


10" 


123877 


35.100 


10" 


334567 


34.767 


10" 


915443 


34.534 


10" 


2520503 


34.210 


10" 


7002043 


33.928 



Let oj{n) represent the number of different prime factors of n. Also, given an integer sequence {wj}^!, 
note that a prime p is said to be a primitive prime factor of if p divides rrii but does not divide any ruj 
for j < i. 

Lemma 2.1.5 (Erdos 1949). Let n be a base 2 pseudoprime. For every k, there exist infinitely many 
squarefree base 2 pseudoprimes with uj{n) — k 

Theorem 2.1.6. There exist infinitely many squarefree base b pseudoprimes n for any b >2 with Lu(n) = k 
distinct prime factors. 
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Proof. Let {nj}°°^^ be an integer sequence of base b pseudoprimes such that each term is greater than its 
preceding term, and U!{ni) — k ~ 1, for any in {nj}^^^. Let pi be one of the primitive prime factors of 
— 1. Since 6"* = 1 (mod pi-Ui) and6'''^^ = 1 (mod p^), -rii is a pseudoprime to base &. We observe 
that b^*^^ = 1 (mod rii) because pi — 1 = (mod (n^ — 1)). As a result, it foUows that = 1 (mod rii). 

Also, = 1 (mod pi ■ Ui) since 6"'^'"^ = • 6"'~^ ■ 6^'"^. Hence, pi ■ rii is squarefree and 

i^{pi ■ ni) = k. Moreover, every integer satisfying pt ■ nt is different because nt is composite, pt > nt, and 
Pi = l (mod {rii - 1)). ■ 

Theorem 2.1.7. For any base b pseudoprime, b > 2, having k > 2 distinct prime factors and for x 
sufficiently large, 

(2.1.8) 3^b,k+i{x)> ^bA^ogf^x). 

Proof. Let n be a pseudoprime with k > \ distinct prime factors. Since n — 1 is the smallest exponent e 
such that p\b'^ — \ and e divides an exponent h such that p | 6'* — 1, it follows from Fermat's little theorem 
that p I IP^^ — 1. Thus, from Zsigmondy's theorem, there exists a prime p > n for which p \ — 1 and 
n — 1 \ p — 1 for b > 2. As a result, 

(2.1.9) np I - 1. 

On the other hand, since np — 1 = n{p — 1) + n — 1 and n— l\p — l,n— l\np— 1 and np \ 6"^"-'^ — 1. If 
we let n, m € N*, the set of positive natural numbers, such that n ^ m and p > n, q > m, then np ^ niq 
for primes p and q. However, suppose we let np = mq and p > n, then m \ p. Hence, m > p and m > n. 
Unfortunately, the latter statement is contradictory, and as a result np ^ mq. If n and m are two different 
base b pseudoprimes with k > 2 distinct prime factors, then np and mq are distinct pseudoprimes as well. 
From ((2T9| . 

(2.1.10) p I (6"^ - 1)(&'^ +1), 
and 

(2.1.11) p<b'^ + l<b^. 

If n < logf^x, then pn < x^log^x < x. It then follows that for every base b pseudoprime n with k distinct 
prime factors, n = piP2 ■ ■ ■ Pk < log^ x, there is at least one base b pseudoprime such that piP2 ■ ■ ■ PkP < x. H 

2.2. Carmichael Numbers. Improving upon Erdos' results in [17], Pomerance [11] sharpened the upper 
bound on C'{x). 

Theorem 2.2.1 (Pomerance 1981). 

(2.2.1) C{x) < X • exp I - ^^1^ I log(3) x + log^^^ 

log^ ' X I 

In the other direction, Alford, Granville, and Pomerance proved a lower bound for C{x) for x sufficiently 
large [15_. 

Theorem 2.2.2 (Alford-Granville-Pomerance 1994). 

(2.2.2) C(x)>a;?, 
thus there are infinitely many Carmichael numbers. 

Recently, Harman improved this lower bound [18] . 
Theorem 2.2.3 (Harman 2005). 

(2.2.3) C(x) > x°-33336704^ 

It is not yet even known if C{x) > x^ . 
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Table 5. Values of /3 



Bound 


10'^ 


10* 




10" 


10^ 


10« 


10« 


10^" 


C{x) 


1 


7 


16 


43 


105 


255 


646 


1547 


P 





0.21127 


0.24082 


0.27224 


0.28874 


0.30082 


0.31225 


0.31895 



Bound 


10'' 


10^^ 


10^^ 


10^4 


10^^ 


10^*^ 


C{x) 


3605 


8241 


19279 


44706 


105212 


246683 


P 


0.32336 


0.32633 


0.32962 


0.33217 


0.33480 


0.33700 



Bound 


10'' 


10^« 


10^'-* 


1020 


10^^ 


C{x) 


585355 


1401644 


3381806 


8220777 


20138200 


P 


0.33926 


0.34148 


0.34364 


0.34575 


0.34781 



We provide in Table [5] a computation of the exponent /3 for which C{x) = x^ for a sufficient value of x 
up to 10^1. 



Conjecture 2.2.4 (Granville-Pomerance 2001). // we let Cslx) be the counting function for Carmichael 
numbers with 3 distinct prime factors, then 



(2.2.4) 

where 
(2.2.5) 

(2.2.6) 
(2.2.7) 



C3{x)^T3 



X3 



T3 



dt 



log'^x 27 J 2 \ogH' 



T3 = K3X « 2100, 



77.1727, 



K3 



E gcd(n, 6) -j-r p 
„4/3 11; 



„4/3 11„_3 
n>l p\n a<b<c, n—abc 

p>3 a, b,c pair wise coprime 



E 5\a,b,c)\{ 



p>3 



{2, ifa = b = c^O (mod 3); 
1, if otherwise. 

and LOa,b.c{p) is the number of distinct residues modulo p represented by a, b, c. 

Recent provisional estimates by Chick and Davies [14] of the slowly converging infinite series K3 suggest 
that K3 = 27.05 which gives = 2087.5. 

3. On the Distribution of Carmichael Numbers 

3.1. Two Conjectures Regarding fc-prime Pseudoprimes and fc-prime Carmichael numbers. We 

conjecture the following relations: 

Conjecture 3.1.1. For any fixed k>2, let ^b.k{x) denote the counting function for base b pseudoprimes 
with k distinct prime factors, and let l3^b{x) denote the counting function for base b pseudoprimes. Asymp- 
totically, 



(3.1.1) 



o(l). 



^bix) 

In other terms, for any fixed base b > 1, the k-prime base b pseudoprimes, ^b,k{x), form a set of relative 
density in the set of all base b pseudoprimes, ^bix), for that same value ofb. 
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We are only able to partially support Conjecture 13. 1.11 First, we express the ratio ^^'^(^^ as. 



(3.1.2) 



i,{x) k{x) 



where the maximum number of distinct prime factors, k{x), of any integer < a; is k{x) ^ i^gf^) x ' ^'^^ ^'^Sb^ 
denote the the j-fold iteration of the base b logarithm. Thus, 

g{x) fc-1 k(x) 

(3.1.3) 5Z'^MW = I]^M(2^) + ^b,fc+ ^bA^)- 

1=2 1=2 i=fc+l 



X 



Due to Theorem[2lL7l for any /i < fc in (|3T3)) . ^b,k{x) > ^6,ft(log[''"''^ x), and for any w > fc in ((3T3)) . 

^bA^) > ^h,fc(logi"'"''^a;). Hence, 

fc-1 

(3.1.4) J2 ^bA^) ^ ^b,2(logl'~'^ x) + ^f,,3(logi'"'^ x) + ■ • • + i3^;,,fc_i(logba;) 

i=2 

We cut off the terms froj 
large to satisfy (jS.l.Sp . 



We cut off the terms from proceeding until ^^°f2)^ because if such were the case, then no x could be sufficiently 



k{x} 

(3.1.5) ^"A^^ ^ ^b,k+i{log,x) + ■■■ + ^6,.(,)(logi'^^"^-'' x), 

i=k+l 

where r{x) is any function that grows slower than log* a;, the iterated logarithm. We explicitly define log* x 
as 



(3.1.6) log*x 



ifa;<l; 

1 + log* (logs) ifa;>l 



Remark 3.1.2. We should note that the support for Conjecture 13 . 1 . II is rather weak. This is largely due to the 
weakness of Szymiczek's construction, <^b.fc+i(x) > =^f,^fc(logj, x), in Theorem l2.1.7l We believe that the latter 
relation can be strengthened if a polynomial decrease can be proven. In other words, if ^f,,fc+i(2;) > ^b,k{x'^) 
for some c £ (0, 1). Similarly, in our proof of Conjecture 13. 1.11 we defined the function r{x) as any function 
that grows slower than log* x, the iterated logarithm. Although it is not hard to see that any function 
growing faster than log* x will fail, it is not obvious whether any function growing at the same rate as log* x 
will succeed. However, we have several reasons to strongly believe that r{x) = log* x. First, for practical 
values of X < 2^^^^^ the iterated logarithm grows much more slowly than the logarithm. Second, the iterated 
logarithm's relation to the super-logarithm also supports its slow growth. Third, higher bases give smaller 
iterated logarithms, and log* x is well defined for any base greater than exp This implies that for any 

base b > 2, the iterated logarithm will grow even more slowly for higher pseudoprime bases. 

Conjecture 3.1.3. For any fixed k > 3, let C'k{x) denote the number of k -prime Carmichael numbers up 
to X, and let C{x) denote the Carmichael number counting function. Asymptotically, 

(3.1.7) |^^ = o(l). 

3.2. Support for Conjecture 13.1.11 and Conjecture 13.1.31 So far, the claim established by Conjec- 
ture [^1111] is not yet borne out by the data in Table [51 We believe that the ratio '^^'^l^^) approach 0, but 

may do so slowly at first. On the other hand, it appears that the ratio in Table [7] rapidly approaches 

0, thereby supporting Conjecture 13. 1.31 
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Table 6. Values of 



Bound 






.0^2 (x) 


10^ 





3 


0.00 


10* 


11 


22 


0.50 


10^ 


34 


78 


0.44 


10*^ 


107 


245 


0.44 


10^ 


311 


750 


0.41 


108 


880 


2057 


0.43 


109 


2455 


5597 


0.44 


lOi" 


6501 


14884 


0.44 


10" 


17207 


38975 


0.44 


10" 


46080 


101629 


0.45 


10" 


123877 


264239 


0.47 


10" 


334567 


687007 


0.49 


10" 


915443 


1801533 


0.51 


10" 


2520503 


4744920 


0.53 


10" 


7002043 


12604009 


0.56 



Table 7. Values of 



Bound 




Cix) 


- ■ , , , 

C{x) 


10^ 


1 


1 


1.00 


10^ 


7 


7 


1.00 


10^ 


12 


16 


0.75 


10^ 


23 


43 


0.53 


10^ 


47 


105 


0.45 


108 


84 


255 


0.33 


109 


172 


646 


0.27 


10" 


335 


1547 


0.22 


10" 


590 


3605 


0.16 


10" 


1000 


8241 


0.12 


10" 


1858 


19279 


0.096 


10" 


3284 


44706 


0.073 


10" 


6083 


105212 


0.058 


10" 


10816 


246683 


0.044 


10" 


19539 


585355 


0.033 


10" 


35586 


1401644 


0.025 


10" 


65309 


3381806 


0.019 


1020 


120625 


8220777 


0.015 


1021 


224763 


20138200 


0.011 



Furthermore, Pomerance, Selfridge, and Wagstaff's famous results [10^ support both conjectures. In 
Conjecture 1 of their paper, they beheve that for each e > 0, there is an xo{e) such that for all x > xo(e), 

-{2 + e}logx-log^^^a;~i 



(3.2.1) C(x) > a;•exp|- 



log^2) 2, J • 

Pomerance, Selfridge, and Wagstaff [TO] show that ^bA^) < Ofc(a;2'=/(2fc+i))^ jf (pXT|) is true, then the 
pseudoprimes "with exactly k prime factors form a set of relative density in the set of all [pseudoprimes]" 
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[TU] , Similarly, in Theorem 7 of Granville and Pomerance [|8j, it is proven that Ck{x) < a;2/3+o,(l)^ g^^^ jf 

(IX^ holds, "then for each k, Cfc(x) = o(C(a;))" ^U\. ' 

Interestingly, we can also support the statements in Conjecture 13.1.11 and Conjecture 13.1.31 by relating 
them to their composite superset. Let the number of composites < x with k distinct prime factors be 
denoted by TTkix) and let the number of composites < x with k prime factors (not necessarily distinct) be 
represented by Tk{x). Hence, we can prove upper and lower bounds for 7rfc(x). In 22.18.2 of Hardy and 
Wright [7J for fc > 1, 



(3.2.2) klTTkix) < Hfc(x) < kWkix), 

where Uk{x) = ^ + 0{j^) in 22.18.5. In 22.18.24, since Mx) = T^kix)\ogx - ^ 
kx{\og^^'^ x)''-'^ for fc > 2 and ^ ^^^dt = 0(x), Hfc(.T) - '"'"'^'"fog^'"-''"' ■ As a resuh, it follows that 

In the same respect, a lower bound for tt/j can be formulated. In 22.18.3 it is proven that, 
(3.2.4) Tk{x) - TTkix) < 1^ J2 l:=Hfe_i(a;). 

Since nkix) > Tk{x) - Ilk-i{x) and -Kkix) > - Hfc_i(a;), 

, , ^ /x(log(!M^\ (fc-lMW^M^ 



dt 



(fc — l)!logx / log a; V^og^; 

We can improve the upper bound given in p.2.3|) to an equality, 

x(log(2^ xf-^ 



(3.2.5) ^k{x) 



(fc — 1)! logx 



By the Erdos-Kac Theorem |21| , we can formulate the probability that a number near x has k distinct prime 
factors using the fact that these numbers are distributed with a mean and variance of log^ ' x. Hence, setting 
log'^' x as the A of the Poisson distribution P(fc; A) and taking its limit for any fixed fc, 

(3.2.6) hm P(fc; A) ^ lim (W^^ ^)-;; exp{- x} ^ 

X — >oo X — >oo (k — 1)! 

where the asymptotic error bound is given by 0( jog(2) ^ ) [13] . However, we caution the reader to consider 
that just because the probability of a general composite near x having fc distinct prime factors goes to 0, 
does not necessarily fully prove that this probability will hold for either 3^b,k{x) or Ck{x). 

3.3. An Alternate Conjecture. From Coniecture l3.1.1l and Coniecture l3.1.3l it is evident that the fc-prime 
pseudoprimes and the fc-prime Carmichael numbers are much more sparsely distributed than the set of all 
pseudoprimes and Carmichael numbers, respectively. We hypothesize that if fc is minimized for both the 
fc-prime pseudoprimes and the fc-prime Carmichael numbers, then the ratios '^^'^i^^-^ and ^jj^ will roughly 
achieve the same values for large enough x. We also recommend using the minimum number of distinct prime 
factors for both the pseudoprimes and the Carmichael numbers because first, there is no overlap between the 
three-prime Carmichael numbers and two-prime pseudoprimes and second, the distinct prime factors cannot 
be arbitrarily chosen. This idea leads us to believe that, 

C(.) . ^3(.)^.(.) 



6.2 



(^) 
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As a result, assuming Coniecture l2.1.31 Coniecture l2.1.41 Coniecture l2.2.41 and Coniecture ll.0.21 the amount 
of Carmichael numbers < x given by the counting function C{x) is conjectured to be for x sufficiently large, 



(3.3.1) 

where 
(3.3.2) 

and 
(3.3.3) 



Cix) 



Ip'xi 



-01X2 log^ X 



dt 



logx ■ L{x) 



L{x) 



^ C 



V''i 



27C' 



In Table m the computed values of i/'' and ip'-^ up to lO^^ are given. Hence, not only does Corollary 11.0.31 

Table 8. Values of -0' and ip[ 



Bound 


C{x) 


69.51x§ 




0' 




log x-L{x) 


L{x) 


10^ 


1 


301.95 


1092.82 


0.2302 


0.0024 




7 


594.43 


2835.17 


0.8185 


0.0063 


10^ 


16 


1316.29 


6640.29 


0.8449 


0.0062 


10^ 


43 


3131.53 


14806.24 


0.9545 


0.0075 


10^ 


105 


7826.17 


32411.27 


0.9326 


0.0083 


108 


255 


20282.91 


71150.56 


0.8739 


0.0092 


109 


646 


54070.80 


159157.24 


0.8305 


0.0104 


lOi" 


1547 


147451.71 


367012.00 


0.7293 


0.0108 


10" 


3605 


409716.38 


878601.38 


0.6116 


0.0105 


1012 


8241 


1156637.85 


2188667.23 


0.4953 


0.0097 


1013 


19279 


3309970.24 


5664006.88 


0.4049 


0.0087 


1014 


44706 


9585268.36 


15162465.67 


0.3242 


0.0076 


1015 


105212 


28049810.91 


41763706.96 


0.2607 


0.0065 


1016 


246683 


82852448.55 


117743387.56 


0.2070 


0.0054 


lOi^ 


585355 


246785788.13 


338238941.70 


0.1649 


0.0044 


1018 


1401644 


740679196.52 


986503770.93 


0.1315 


0.0037 


1019 


3381806 


2238429061.23 


2913197684.15 


0.1050 


0.0030 


1020 


8220777 


6807841639.58 


8692508977.60 


0.0839 


0.0024 


1021 


20138200 


20826296835.28 


26167265004.43 


0.0672 


0.0020 



fit the proven bounds for C{x) given in Theorem 12.2.11 and Theorem 12.2.31 but both ip' and -05^ appear to 
be approaching constant values. However, there are several reasons as to why Corollarv 11.0.31 mav not be 
necessarily borne out by the data in the above table. For instance, the infinite series K3 is slowly convergent, 
and it is not until 102^ that K3 appears to approach its estimated value of 2087.5. However, the primary source 
of inaccuracy is due to Conjecture 12. 1 .31 Since Pomerance's conjecture for the distribution of pseudoprimes 
is applicable for sufficiently large x and pseudoprime counts have only recently been conducted to lOi'' by 
Galway and Feitsma, we are not sure how "sufficiently large" x must be for Coniecture l2.1.3l to be an accurate 
model for pseudoprime distribution. Lastly, x must also be immensely large in order for '^^^^^^•^ = o(l). 

3.4. An Improved Heuristic Argument. As mentioned before, Pomerance's heuristic arguments sup- 
porting his conjecture in (|2.2.2p involve the distribution of smooth numbers. And, if secondary terms 
exist, then it would be worthwhile to sharpen these heuristics to produce a viable conjecture for C{x). 
Let ^!{x,y) denote the number of y-smooth numbers < x and let '^'{x,y) denote the number of primes 
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p < X for which p — 1 is squarefree and its prime factors are < y [lOj . It is conjectured in [llj that for 
exp{i(loga;)^} < y < exp{(loga;)5}, 

(3.4.1) i^r(a,,y)^ 1 vl/'(a;,y). 

X Tr[x) 

If < a < 1, it is well-known [5^ that 

(3.4.2) ^ (a;,exp{c(loga;)"(log(^^ xf}'^ = a:cxp{-{(l -a)/c + o{l)}{\ogx)^-°'{\og'^^^ x)^-''^}. 
Concerning Carmichael numbers, we are interested in the case for which a ~ ^, /3 — 0, and c = 1. Hence, 

(3.4.3) 4' (a;,exp{(logx)5}^ = a;exp{-{l/2 + o(l)}(logx)^(log(^> x)}. 
From (|3.4.ip and (|3.4.3p . we make the following 

Conjecture 3.4.1. For exp{i(loga;)^} < y < exp{(loga;)^}, 

(3.4.4) "^'(x^y) = 7r(a;)exp{-{l/2 + o(l)}(loga;)5(log(2)x)}. 

Let A{x) denote the product of the primes p < loga;/(log''^^ a;)^. Thus, A{x) < x^^^°^'^^^^ as in [10]. If 
we allow ri , . . . , to be the primes in the interval (loga;/(log(^) x)^, (logx)'°s<" ^^^'j with r, - 1 | A{x). By 



Conjecture 13.4. II we have for a; sufficiently large, 

(3.4.5) q = Tr ((log a;)^°s'" exp{-{l + o(l)} log^^) x log^^^ a;}. 

Let mi, . . . , niN be the squarefree composite integers < x composed of ri and let I = loga;/(log'^-' a;)" 



As discussed in [10], we have 

(3.4.6) > ( ] ) > 



As a result, 



(3.4.7) N > 



' (W^) x)^7r((loga;)'°s'" ") exp{-{l + o(l)} log^^^ a; ^ x iogx/(iog< ' x) 

logx j 



Since Euler's function and Carmichael's A function are virtually the same, the lower bound in p.4.7p should 
be applicable to C{x). In fact, from the values of a{x) in Tableland the precision of Conjecture 13.4.11 we 
have reason to believe that this result is asymptotically close to the actual value of C (x) . 

4. One- Parameter Quadratic-Base Pseudoprimes: A Sidenote 

As mentioned earlier, the discovery of Carmichael numbers demonstrated the fallability of Fermat's pri- 
mality test therefore lending to the development of efficient probabilistic primality tests. Baillie, Pomerance, 
Selfridge, and Wagstaff [10] [25] have determined a primality test that is an amalgamation of the Miller-Rabin 
test and a Lucas test. However, even though Pomerance j26j presented a heuristic argument that the number 
of counter-examples up to x was ^ x^~^ for e > 0, we have not been able to find any counter-examples up 
to 10^^. In fact, no precise probability of error has been given about this test either 24J. 

Grantham 27J has also provided a probable prime test known as the RQFT that has a known worst-case 
probability of error of 1/7710 per iteration. 

An even stronger test known as the One-Parameter Quadratic-Base Test (OPQBT) has been given by 
Zhang [24J, and is a version of the Baillie-PSW test that not only has a known probability of error but is more 
efficient than the RQFT except for a thin set of cases. We let u(^ ±2) e Z, let r„ = T (mod - uT + 1), 
and define the ring associated with parameter u as 

Ru Z[T]/{T^ - uT + 1) = {a + bTu : a,b e Z}. 
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We then define an odd integer n > 1 as an OPQBT pseudoprime for < u < n with 
where in the ring n must pass 

(4.0.8) T^-' = 1 (mod n). 

Moreover, n is defined as an OPQBT strong pseudoprime if for some i = 0, 1, • ■ • ,k — 1, either 
(4.0.9) Tg" = 1 (mod n), 

or 

(4.0.10) « = -1 (mod n), 

in which for q odd, n~ e — 2^q 

We have verified that there are no OPQBT pseudoprimes up to 10^^. Let the counting function ff{x) 
denote that number of OPQBT pseudoprimes < x and let ,S^ff{x) denote the number of strong OPQBT 
pseudoprimes < x. The best upper bound we are able to prove is 

(4.0.11) yff{x) < e(x) < X ■ L{x)^, 

since an upper bound on the pseudoprimes is applicable to an upper bound on the OPQBT pseudoprimes 
and strong OPQBT pseudoprimes. 

Based upon Erdos' construction [17] and Pomerance's heuristics [26], in the interval , for any fixed 

j > 4 and H sufficiently large, there are most likely exp{iJ-^(l — 4/j)} counter-examples to Zhang's primality 
test, meaning that there are at least x^~'^^^ counter-examples below x = expjff^}. Thus, for arbitrary j, 
the number of counter-examples to the OPQBT becomes generalized to ^ x^~'^ for e > 0. In other words, 
there are infinitely many counter-examples to Zhang's OPQBT. 

5. Acknowledgements 

I would like to express my gratitude to Charles R. Greathouse IV who provided invaluable guidance in 
the direction of this paper; Professor Carl Pomerance who offered helpful comments on this paper in its 
initial stages; Harvey Dubner who provided guidance on the experimental aspects of this paper; Johan B. 
Henkens who helped me count the base 2 pseudoprimes and 2-strong pseudoprimes up to lO^'' using Galway 
and Feitsma's data; Professor Kazimierz Szymiczek who clarified Theorem 12.1.71 and David H. Low who 
assisted me with formatting and typesetting errors. 1 would also like to thank Professor William F. Galway 
for sharing his viewpoints regarding the pseudoprimes with k distinct prime factors. 



ON THE DISTRIBUTION OF CARMICHAEL NUMBERS 



15 



References 

1. C. Pomerance and R. Crandall, Prime Numbers: A Computational Perspective, New York: Springer, 2005. 

2. P. Erdos, On the Converse of Fermat's Theorem, Amer. Math. Monthly, 56 (1949): 623-24. 

3. W. F. Galway, Research Statement, 2004, http: / /www. math. uiuc.edu/~galway/resear ch-statement.pdf| 

4. W. F. Galway, Tables of pseudoprimes and related data, 2002, |http: / /oldweb.cecm.sfu.ca/pseudoprime/psplelS.gz 

5. W. F. Galway, The Pseudoprimes below 2''*, Simon Fraser University, 2002, 
jhttp: / /oldweb.cecm.sfu.ca/pscudoprime/psp-search-slidos.pdf 

6. C. Pomerance and D. M. Gordon, The Distribution of Lucas and Elliptic Pseudoprimes, Math. Comp. 57 (1991): 825-38. 

7. E. M. Wright and G. H. Hardy, An Introduction to the Theory of Numbers, Oxford: Clarendon P, Oxford UP, 1998. 

8. A. Granville and C. Pomerance, Two Contradictory Conjectures Concerning Carmichael Numbers, Math. Comp. 71 (2001): 
883-908. 

9. C. Pomerance, A New Lower Bound for the Pseudoprime Counting Function, Illinois J. Math. 26 (1982): 4—9. 

10. C. Pomerance, J. L. Wagstaff, and S. S. Wagstaff, jr, Pseudoprimes to 25 • 10^, Math. Comp. 35 (1980): 1003-026. 

11. C. Pomerance, On the Distribution of Pseudoprimes, Math. Comp. 37 (1981): 587—93. 

12. K. Szymiczek, On Pseudoprimes which are Products of Distinct Primes, Amer. Math. Monthly 74 (1967): 35—37. 

13. G. Tenenbaum, Introduction to Analytic and Probabilistic Number Theory, Cambridge: Cambridge UP, 1995. 

14. G. H. Davies and J. M. Chick, The Evaluation of K3, Math. Comp. 77 (2008): 547-550. 

15. A. Granville, C. Pomerance, and W. R. Alford, There are Infinitely Many Carmichael Numbers, Ann. of Math. 140 (1994): 
703-722. 

16. R. G.E. Pinch, The Carmichael Numbers up to 10^^, Proceedings Conference on Algorithmic Number Theory, Turku, May 
2007. Turku Centre for Computer Science General Publications 46, edited by Anne-Maria Ernvall-Hytonen, Matti Jutila, 
Juhani Karhumaki and Arto Lepisto. 

17. P. Erdos, On pseudoprimes and Carmichael numbers, Publ. Math. Debrecen 4 (1956): 201—206. 

18. G. Harman, On the number of Carmichael numbers up to x. Bull. Lond. Math. Soc. 37 (2005): 641-650. 

19. R. G.E. Pinch, The Carmichael Numbers up to 10 to the 21, Eighth Algorithmic Number Theory Symposium ANTS- VIII 
May 17-22, 2008 Banff Centre, Banff, Alberta (Canada). 

20. A. Korselt, Probleme chinois, L'intermediaire de mathematiciens 6 (1899): 142—143. 

21. P. Erdos and M. Kac, The Gaussian law of errors in the theory of additive number theoretic functions, Amer. J. Math 62 
(1940): 738-742. 

22. R. D. Carmichael, Note on a new number theory function. Bull. Amer. Math. Soc. 16 (1910): 232-238. 

23. D. Shanks, Solved and unsolved problems in number theory, 3rd ed., Chelsea, New York, 1985. 

24. Z. Zhang, A one-parameter quadratic- base version of the Baillie-PSW probable prime test. Math. Comp. 71 (2002): 1699- 
1734. 

25. R. Baillie and S. S. Wagstaff, jr., Lucas pseudoprimes. Math. Comp. 35 (1980): 1391-1417. 

26. C. Pomerance, Are there counter-examples to the Baillie-PSW primality test?, Dopo Le Parole aangeboden aan Dr. A. K. 
Lenstra (H. W. Lenstra, jr., J. K. Lenstraand, P. Van Emde Boas, eds.), Amsterdam, 1984. 

27. J. Grantham, A probable prime test with high confidence, J. Number Theory, 72 (1998) 32-47. 

28. Emails exchanged between A. Nayebi and C. Pomerance. 

29. A. Granville, Smooth numbers: computational number theory and beyond. Mathematical Sciences Research Institute Pub- 
lications, 44 (2008) 1-58. |http:// www. math. leidenuniv.nl/l psh / ANTproc / 09andrew.pdf 

E-mail address: aran.nayebiOgmail . com 



